Electron Glass in Ultrathin Granular Al Films at Low Temperatures 
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Quench-condensed granular Al films, with normal-state sheet resistance close to 10 kfi/O, display 
strong hysteresis and ultraslow, non-exponential relaxation in the resistance when temperature is 
varied below 300 mK. The hysteresis is nonlinear and can be suppressed by a dc bias voltage. The 
relaxation time does not obey the Arrhenius form, indicating the existence of a broad distribution of 
low energy barriers. Furthermore, large resistance fluctuations, having a l//-type power spectrum 
with a low-frequency cut-off, are observed at low temperatures. With decreasing temperature, 
the amplitude of the fluctuation increases and the cut-off frequency decreases. These observations 
combine to provide a coherent picture that there exists a new glassy electron state in ultrathin 
granular Al films, with a growing correlation length at low temperatures. 



PACS numbers: 72.80.Ng, 73.50.-h, 64.60. My, 05.40, 



The electron glass was first predicted to exist in dis- 
ordered interacting systems nearly two decades ago [[!). 
Such nonergodic behavior is very interesting because one 
normally expects electron systems to relax rather rapidly. 
Over the years, a number of studies have been reported 
in which electrons display glassy dynamics that are of- 
ten associated with non-exponential relaxation extend- 
ing over many decades in time, such as the field-effect 
conductance measurements in compensated GaAs 
amorphous indium-oxide films ||, and ultrathin Bi/Gc 
and Pb/Ge films [jjj. The glassy behavior is believed 
to arise from the electron-electron (e-e) interactions and 
the Coulomb gap ^-0. Recently, the electron glass has 
received renewed interest || as the subject of e-e inter- 
actions has become a central topic in understanding the 
metal-insulator transition in two-dimensions Q. How- 
ever, the precise role that the Coulomb gap plays in the 
observed glassy behavior is not clear since there is no si- 
multaneous measurement of the conductance relaxation 
and the single particle density of states. In addition, 
there is no direct experimental effort to probe the corre- 
lation length in the glassy phase. 

In this Letter we report glassy behavior in the normal 
state of quench-condensed weakly insulating granular Al 
films of sheet resistance, Rn, of about 10 kfi/D at 300 
mK. We measure the relaxation of Rn after the temper- 
ature was varied. We have focused on weakly insulat- 
ing films because Rn on the order of 10 kfl/D is easy 
to measure using sensitive ac lock-in techniques, which 
turn out to be crucial in measuring the resistance fluc- 
tuations described below. We observed that, below 300 
mK, the resistance was strongly hysteretic and displayed 
ultraslow, non-exponential relaxation as the temperature 
was varied. We have also observed strong nonlinear be- 
havior in the hysteretic regime. What was unique to our 
work was the first observation in quench-condensed metal 
films of large resistance fluctuations below 100 mK. We 



argue that these observations indicate the existence of a 
glassy electron state in ultrathin granular Al films with 
a growing correlation length at low temperatures. 

Our Al films were quench-condensed onto glass sub- 
strates using 99.999% purity Al sources in UHV condi- 
tions inside a dilution refrigerator, with the substrates 
being held near 20 K during evaporation. The film thick- 
ness was near 25 A. After warming up to room temper- 
ature, the films showed a granular morphology with a 
typical grain size of about 300 A, as seen m scanning 
force microscopy studies. The films had a multi-lead pat- 
tern with an area of 3x3 mm 2 between the neighboring 
leads. An analog lock-in amplifier, operating at 27 Hz 
with a time constant of 3 seconds, was used to measure 
the four-probe ac resistance. The ac probe current was 
fixed at 1 nA, producing an ac voltage of about 10 fj,V 
across a film section. For nonlinear studies described be- 
low, an additional dc bias voltage was applied. All the 
data described below were measured in the normal state, 
with superconductivity being suppressed by a magnetic 
field above the spin-paramag netic limit f(|0 of about 
48 kG. To date, we have performed detailed studies on 
films of Rn ~ 10 kf2/D. Such films appeared to be very 
uniform, with R a varying less than 5% among the var- 
ious sections of the multi-lead pattern. These films are 
far above the percolation threshold since Rn scales with 
inverse film thickness. 

In the inset to Fig. 1, we show a typical resistance 
hysteresis measured on one Al film as temperature was 
cycled between 500 and 30 mK, in the high-field nor- 
mal state with a field, H± = 52 kG, applied perpendic- 
ular to the film plane. Throughout all our experiments, 
the rates of cooling and heating were kept constant on a 
logarithmic temperature scale, and the cooling and the 
heating cycles each took 1.5 hr. We can exclude tem- 
perature lag or heating as the cause of the hysteresis, 
because we have not observed any hystersis in Li films 
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FIG. 1. Main figure shows time constant r vs. 1/T for 
various applied fields: perpendicular 52 kG (up triangles), 
parallel 52 kG (down triangles), perpendicular 90 kG (cir- 
cles), and parallel 90 kG (diamonds). Inset shows a typical 
resistance hysteresis, with the arrows indicating the direction 
of the temperature sweeps. 



of Ro ~ 10 ktt/a and Be film of Rn = 3 ~ 500 kfl/U, 
both quench-condensed on the same glass substrates. For 
relaxation studies, the films were first cooled from 500 
mK to the desired temperature, T, using the above log- 
arithmic rate, at which R a was measured as a function 
of time. In Fig. 2, we show the Rn vs. time curves 
for a number of fixed T in a field of H± = 52 kG. We 
found that the curves at higher T, such as 100 mK, 
could be fit very well to a stretched-exponential form, 
[Ra(cx>) - Ra(t)]/[Ra(oo) - Rn(0)] = exp[-(t/T)"<}, over 
three decades in time. With decreasing T, the time con- 
stant r increased sharply. The exponent 7 scattered be- 
tween 0.6 and 0.8 at low T without a clear trend. We note 
that below 100 mK, the large resistance fluctuations seen 
in Fig. 2 made it difficult to accurately determine r and 
7. In the main part of Fig. 1, we plot r vs. 1/T for two 
field values in both perpendicular and parallel field orien- 
tations. Data in Fig.l appear to fall into two groups: one 
at 52 kG and the other at 90 kG, suggesting that r was 
larger in higher fields and was insensitive to field orien- 
tation. The almost linear dependence of r on 1/T in Fig. 
1 indicates that, with decreasing T, r increased much 
slower than the Arrhenius law, t(T) = TQexp(E a /T), 
where E a is known as the activation energy describing 
the typical energy barrier to relaxation [ |12| . Such non- 
Arrhenius behavior can be explained only if there exists 
a very broad distribution of low-energy barriers. Relax- 
ation occurs over lower energy barriers with decreasing 
T, leading to a relaxation time that increases slower than 
the Arrhenius law. 

The above conclusion is also supported by the nonlin- 
ear effects in the hysteresis. For nonlinear studies, we 
have measured Rn in the presence of a dc bias voltage, 
Vbias, across the films, and have observed that the hys- 
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FIG. 2. Time traces of film sheet resistance measured at 
temperatures as labeled on the graph, in a field of H± — 52 
kG. Smooth lines are fits to a stretched-exponential form. 



teresis can be suppressed if Vbias is large enough. The 
experimentally relevant Vbias was of the order of 0.5 mV. 
It leads to a bias current of 0.05 /iA and a Joule heating of 
2.5 xl0~ n W. We believe that the nonlinear effects were 
not caused by heating. Consider the hysteresis shown in 
the inset to Fig. 1, we observed that, between 100 mK 
and 200 mK, the heating curve shifted downward and the 
cooling curve shifted upward with increasing Vbias- It is 
obvious that had heating by the bias current been signifi- 
cant the cooling curve should have also shifted downward 
with increasing Vbias, since the resistance of weakly in- 
sulating films should decrease with increasing T. In Fig. 
3(a), we plot the width of the hysteresis, AR, versus T 
for various Vu as values. Such data are obtained by sub- 
tracting the cooling curve from the heating curve in a 
thermal cycle such as the one shown in the inset to Fig. 
1. It is quite revealing in that the bias voltage does not 
uniformly suppress the hysteresis across the entire range 
of T. Instead, a small Vbias suppresses the low-T com- 
ponents of the hysteresis only. With increasing Vbi as , 
progressively higher-T components of the hysteresis are 
suppressed. We believe that relaxation processes over 
low-energy barriers are suppressed by the bias voltage 
in the same way the low-T components of the hystere- 
sis are suppressed. Nevertheless, we have yet to develop 
quantitative methods to analyze the data in Fig. 3(a) to 
obtain the distribution of low-energy barriers. We note 
that the falling off of AR at low T in the Vbias = curve 
in Fig. 3(a) was an experimental artifact because we did 
not wait at 30 mK for the system to relax after cooling 
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FIG. 3. (a) shows the hysteresis width as a function of tem- 
perature for H± = 52 kG. Curves from top to bottom corre- 
spond to Vtias = 0, 0.3, 0.4, 0.5, 0.98, 1.5 3.5, and 10 mV, 
respectively, (b) shows the scaled hysteresis width vs. Vbias 
at 175 mK. There are four types of symbols, corresponding 
to applied fields of 52 kG or 90 kG, in either perpendicular or 
parallel orientations, (c) shows a log-log plot of the tempera- 
ture dependence of the threshold voltage. The solid line is a 
power- law fit with a power of 3.1 ± 0.3. 



down from 500 mK. Had we waited long enough before 
heating up towards 500 mK, AR would have been the 
largest at 30 mK for the Vbias = curve. 

If we make a vertical cut at a fixed T, such as 175 mK, 
in Fig. 3(a), we obtain AR as a function of Vu a s- We 
have observed that the absolute value of AR varied with 
the strength and the orientation of the magnetic field. 
However, if we scale AR by its value in the Vbias = limit, 
Ai?(0), the resulting scaled width, AR/AR(0), appears 
to be independent of the strength and the orientation of 
the field. This is demonstrated by the good data collapse 
in Fig. 3(b) in which we plot four sets of such scaled 
data at 175 mK as a function of Vbias for two field values 
in both perpendicular and parallel orientations. In this 
semi-log plot, AR falls off steeply for Vu a s > 0.7 mV. 
This suggests that there is a field-independent threshold 
bias, Vq ~ 0.7 mV, that suppresses the hysteresis. By 
making vertical cuts at various temperatures in Fig. 3(a), 
we observe that Vq decreases with decreasing T. This is 
consistent with the picture that relaxation occurs over 
lower energy barriers at lower T; as a result, the value 
of Vq for suppressing the hysteresis is also lowered. We 
plot Vq vs. T in Fig. 3(c). The flattening of the data 
below 100 mK is again due to the experimental artifact 
discussed earlier. For the high-T part of the data, Vq 
can be fit to a power law temperature dependence, Vq ~ 
T 31±0 - 3 , as shown by the solid line in the inset to Fig. 3. 

The results discussed above suggest the existence of a 



glassy state in our quench-condensed ultrathin granular 
Al films. The low temperature scale of our experiments 
excludes the possibility of glassy structural relaxation, 
such as atomic diffusion between metastable configura- 
tions similar to those observed in a-Si:H |13| . Tunnel- 
ing studies have shown in such Al films Jl4|, as well as 
in other films of similar sheet resistance Jl5|, that an 
anomaly in the density of state exists at the Fermi en- 
ergy. However, a true Coulomb gap has not been found 
in such weakly insulating films. Thus it is not clear what 
role, if any, the Coulomb gap plays in the observed glassy 
relaxation. In addition, we did not find any hysteresis in 
quench-condensed Li films of R n ~ 10 kSl / □ and Be films 
of R n = 3 ~ 500 kfi/D, even though a true Coulomb gap 
has been observed in Be films |fl6|| . While we were not 
able to investigate the morphology of the Li films, which 
become unstable in air, scanning force microscopy stud- 
ies of the Be films warmed up to room temperature did 
not find any granular structure down to 1 nm. Because 
of the morphological difference, we propose that granu- 
larity plays an important role in the glassy behavior of 
our Al films. 

Consider a simple conduction model for granular met- 
als in which an electron moves from one neutral grain 
to another nearby one to create a charge- ant icharge pair 
Jr7[ . The energy cost for creating such a pair is 2E C , 
where E c = e 2 /2C is the grain charging energy, with C 
being the capacitance of a grain. E c can be significant 
if the grains are small, and can vary strongly with loca- 
tion due to the randomness in grain sizes. This results 
in a rough potential background for the charge carriers. 
Transport is mediated via the ionization of such pairs 
p8| , with the energy associated with ionization being 
the Coulomb attraction between the pair, which, in two 
dimensions, has a logarithmic form with a cutoff length 
|fl9f . This model has motivated experiments searching for 
a finitc-T Kosterlitz-Thouless-Berezinskii charge unbind- 
ing transition in arrays of Josephson junctions |l9|^o[] and 
granular films |2lf| . However, in a recent theoretical study 
of capacitively coupled grains, Granato and Kosterlitz 
p2| have found that finite-T transitions are suppressed 
by disorder. They found instead a charge glass with a 
correlation length that diverges as temperature decreases 
to zero. They have also predicted that nonlinear behav- 
ior sets in at a characteristic voltage V c ~ T 1+n , with a 
thermal critical exponent n ks 1.7 characterizing the T 
= charge glass transition. It is interesting to note the 
closeness of this predicted power-law dependence with 
the power-law fit in Fig. 3(c). However, we caution that 
our fit has a limited range. Within this picture of granu- 
lar transport, we can qualitatively understand the behav- 
ior of the films in an applied magnetic field as shown in 
the main part of Fig. 1. The field splits the singlet states, 
effectively reducing the density of states. As a result, the 
relaxation time increases with increasing magnetic field, 
however, it is insensitive to the orientation of the field. 
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FIG. 4. Power spectra of the resistance fluctuations are 
shown for 30, 60, 80, 100, and 120 mK, corresponding to 
curves from top to bottom, respectively. Four of the five 
curves are shifted by factors labeled next to the curves. The 
straight line shows the l//-type frequency dependence. 

We note that the data in Fig. 2 show drastic increases 
in the size and decreases in the characteristic frequency 
of the resistance fluctuations with decreasing T, as shown 
by the 60 mK and 30 mK curves. This is a unique fea- 
ture in our data that has not been reported in similar 
quench-condensed films. We have analyzed such noise 
data using a fast Fourier transform (FFT). We first ob- 
tain the noise time traces by subtracting the stretchcd- 
exponential fitting curves from the Ru vs. time traces. 
The noise time traces are then analyzed using FFT. Fig- 
ure 4 shows a set of such spectra obtained from data at 
a number of temperatures with H±_ = 52 kG, including 
the three traces shown in Fig. 2. There are two distin- 
guishing features in the spectra. The first is that, the 
T = 30 mK spectrum, which has the broadest frequency 
range, clearly has a 1/ /-type form [^3| over a wide range 
of frequency, indicating that fluctuations occur at all fre- 
quencies. The second is that all the spectra in Fig. 4 level 
off at low frequencies, indicating a low-frequency cut-off 
in the spectra which clearly shifts to lower frequencies 
as T is lowered. The increase in the magnitude of the 
low-frequency fluctuations and the decrease in the cut-off 
frequency with decreasing T in Fig. 4 indicate that the 
number of low-energy states increases. We suggest that 
these two features arise from the collective hopping of 
many correlated electrons. The cut-off frequency should 
be a measure of the characteristic frequency of fluctu- 
ations occurring on the length scale of the correlation 
length. A decreasing cut-off frequency with decreasing T 
indicates an increasing correlation length. We also point 
out that the power spectra for all the temperatures shown 
in Fig. 4 nearly follow a single curve for / > 10 -2 Hz, 
with a frequency dependence of i/j 1 - 7 ^ 08 . Such be- 
havior deviates significantly from the 1// behavior seen 
in the 30 mK curve at low frequencies. Although the ori- 



gin of such behavior is unknown, it is possible that this 
nearly temperature independent high-frequency feature 
is not related to collective electron hopping. Rather, it 
is the emerging 1// part of the spectra at low tempera- 
tures that represents the appearance of collective electron 
hopping. In general, we expect the following relation be- 
tween the cut-off frequency, f c , the correlation length, 
£, and T: f c ~ £~ z ~ T vz , where v is the correlation 
length exponent and z is the dynamical critical exponent. 
Ideally, one measures the noise using a sensitive bridge 
setup in which the large and non-fluctuating resis- 
tance background is eliminated, rather than the simple 
four-terminal methods used in this study. In addition, 
the fluctuations should be measured when the system is 
in a steady-state, instead of in the initial relaxation pro- 
cess as shown in Fig. 2. Future experiments taking the 
above mentioned concerns into account should provide an 
accurate determination of f c and the critical exponents. 

In conclusion, we have found strong evidence for a new 
glassy electron state in quench-condensed ultrathin gran- 
ular Al films, with a growing correlation length at low 
temperatures. We gratefully acknowledge numerous in- 
valuable discussions with S. Teitel, Y. Shapir, Y. Gao, 
and P. Adams. We thank S. Zorba and Y. Gao who per- 
formed scanning force microscopy studies of the quench- 
condensed Al and Be films. 
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